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ABSTRACT 

v  The  variation  in  range  between  two  satellites  in  different  circular 
orbits  is  of  interest  in  many  space  system  problems.  This  report  provides  a 
simple  expression  for  the  range  and  presents  an  analysis  of  the  frequency  and 
duration  of  close  approaches.  The  results  are  applied  to  a  simple  satellite 
search  problem.  Some  relationships  for  the  range  rate  are  also  given.— 
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I.  INTRODUCTION 

The  purpose  of  this  report  is  to  describe  the  behavior  of  the  range 
between  two  spacecraft  in  different  circular  orbits.  Specifically,  we  have  a 
spacecraft,  denoted  by  subscript  0,  in  a  circular  orbit  of  radius  and 
corresponding  period  T^,  which  carries  a  sensor  or  other  system  characterized 
by  an  effective  radius  Rg .  What  is  the  maximum  amount  of  time  which  can  take 
place  before  a  spacecraft,  denoted  by  subscript  1,  with  orbital  radius  R^, 
orbital  period  T^,  and  relative  inclination  I,  will  come  within  the  sphere 
defined  by  Rg?  What  is  the  average  amount  of  time  that  the  spacecraft  will 
spend  within  the  range  Rg?  What  is  the  maximum  rate-of-change  of  the  range 

between  the  spacecraft? 

In  Section  II  a  convenient  form  for  the  range  between  the  two  spacecraft 
is  derived.  Section  III  defines  a  close  approach,  and  gives  the  time 
interval  between  two  such  approaches.  We  find  a  characteristic  interval 
within  which  a  close  approach  will  occur,  irrespective  of  the  relative 
inclination  between  the  two  orbit  planes.  Section  IV  presents  a  calculation 
of  the  average  time  duration  that  the  two  spacecraft  stay  within  a  specified 
range,  and  it  is  observed  that  this  time  is  minimized  when  the  two  orbit 
planes  are  perpendicular.  Section  V  describes  an  application  of  the  results 
to  a  spaceborne  search  problem.  Finally,  Section  VI  presents  a  simple  upper 
bound  on  the  maximum  range  rate  between  the  two  spacecraft. 


II.  RANGE  BETWEEN  SPACECRAFT 


Figure  1  illustrates  the  geometry  of  interest.  The  plane  of  spacecraft  1 
is  inclined  with  respect  to  the  plane  of  spacecraft  0  by  relative  inclination 
angle  I.  This  inclination  angle  is  uniquely  defined  by  requiring  that  for 
I  ■  0,  both  spacecraft  are  moving  in  the  same  direction. 

This  section  presents  a  derivation  for  an  expression  giving  the  distance 
between  the  two  spacecraft.  Consider  two  rectangular  coordinate  systems  with 
the  x-axis  in  each  corresponding  to  the  line  of  intersection  between  the  planes 
of  the  two  orbits  as  shown  in  Fig.  1.  The  coordinate  system  0  will  be  in  the 
plane  of  motion  of  satellite  0  while  coordinate  system  1  will  be  in  the  plane 
of  motion  of  satellite  system  1.  In  each  of  the  planes,  the  corresponding 
equations  of  motion  of  the  satellites  are 


xo  = 


Ro  sinV 


Yo  = 


Rq  COSO)gt 


(II-l) 


X1  *  R1  sin(CJlt  +  ♦) 

Y1  =  R1  cos^it  +  ♦)  (II-2) 

where 


0) 


i 


(II-3) 


is  the  equivalent  angular  rate  of  the  i satellite  in  its  plane  of  motion. 

We  have  chosen  t  *  0  to  be  the  time  that  satellite  0  is  at  the  intersection  of 
the  planes.  is  the  angular  distance  from  satellite  1  to  the  intersection  of 


II-l 


the  planes,  in  plane  1,  at  t  *  0.  Without  loss  of  generality,  we  assume 

T1  >  V 

Now  let  Xj^,  Y^",  Z  '  be  the  coordinates  of  satellite  1  expressed  in  the 
plane  of  satellite  0.  Then: 


X1  “X1 


*  Y^  cos  I 


Zx  =  sin  I 


(II-4) 


Denoting  the  distance  between  the  two  spacecraft  as  d: 


d2  -  (Xq  -  X/)2  +  (Yq  -  Y1")2  +  (0  -  Z/)2 


(H-5) 


Substituting,  and  using 

1  +  cosX  =*  2cos2  j  , 
1  -  cosX  =  2sin2  , 


(II-6) 


(H-7) 


we  obtain 


d2  -  (Rj^-Rq)2  +  4RqR. 


2  I 

1  COS  2 


I  n,  2  /  V-V-M  ,  ,  2  I  2  J0t4wl^ 

J  sln  l 2 /  +  sin  1  cos  \ 2 


(II-8) 


As  a  simple  check,  Eq.  (II-8)  shows  that  the  range  varies  between  Rq  +  R^ 
and  Rq  -  R^  when  I  is  either  0  or  n. 


III.  CLOSE  APPROACHES 


For  convenience,  introduce: 


the  difference  frequency  u)^  A 


(III-l) 


the  sum  frequency 


A  O)0  +  ^ 


(HI-2) 


the  offset  time 


t  A  t  -  — 

=  0)A 


(IH-3) 


In  these  terms  Eq.  (II-8)  can  be  rewritten 

2  2  I  21  2  ^A1 

d  =  (R^“Rq)  +  4R0R1  cos  2  sin  ~Y~ 


2  1  2 
+  sin  cos 


where  the  equivalent  angular  advance  f|  is: 


(III-4) 


n  =  4>(i  +  -±) 

^A 


(HI-5) 


In  Eq.  (III-4),  the  [  ]  contains  a  term  with  the  difference  frequency  period 


T*  = 


Vl 


A  Tl-T0 


(III-6) 


corresponding  to  the  difference  frequency,  and  a  term  with  a  sum  frequency 
period 


Vi 

Ti+To 


(III-7) 


III-l 


k  rri 


*  2 

corresponding  to  Qualitatively,  d  consists  of  a  constant  term,  a  slowly 

varying  positive  term  with  period  T^,  and  a  rapidly  varying  positive  term  with 
period  T^.  (Note  that  is  greater  than  Tq  or  ,  while  is  less  than  either 

To  or  Tr> 

Since  both  terms  within  [  ]  in  Eq.  (III-4)  are  positive,  it  is  clear  that 
the  closest  possible  approach  (i.e.,  d  -  ^”^0^  occurs  on^Y  for  those  T  such 

2£+l 

that  =  kTT  for  k  integral,  and  u^T  +  n  =  (-^ — )7T  for  £  integral.  Depending 
on  the  particular  u)^  and  r\9  the  time  between  such  closest  possible  approa¬ 
ches  could  be  extremely  long. 

However,  as  a  practical  matter,  because  of  the  great  difference  between 
and  0)^  for  typical  cases,  we  will  consider  the  basic  time  interval  T^. 

During  each  (half-open)  time  interval  of  duration  T.,  the  term  sin^  in 

2 

Eq.  (III-4)  has  one  (and  only  one)  zero.  During  the  same  interval,  the  term 

2  +  n  **  2 

cos  - 2 -  will  have  uu/o).  zeros.  Thus  there  will  be  a  zero  of  the  cos 

1  L  A  2 

term  occurring  within  T ^  of  the  time  at  which  the  sin  term  is  zero.  There¬ 

fore,  within  every  time  interval  of  duration  T^,  the  satellites  will  come 
within  the  distance: 


dmin  <  <VR0)2  +  4R0R1  COs2  I  sl°2  1  ^ 


(III-8) 


One  implication  of  this  expression  is  that  the  relative  inclination  between 
the  orbital  planes  has  a  very  small  effect  upon  the  range  at  closest  approach 
for  small  (Small  CJ^/a)j.  corresponds  to  the  case  where  the  satellite 

periods  are  close  to  one  another.)  Furthermore,  the  characteristic  period 
V  within  which  a  close  approach  is  certain  to  occur,  is  unaffected  by  the 
relative  inclination. 


*  2  2 

From  (II-6)  and  (II-7),  the  sin  and  cos  terms  have  equivalent  frequencies 

equal  to  twice  their  arguments. 

**  2 

For  example,  with  10%  difference  in  periods,  the  cos  term  will  have  21  zeros 

in  this  interval. 


IV,  DURATION  OF  CLOSE  APPROACHES 


In  this  section  the  fraction  of  the  basic  time  interval  T.  that  the  dis- 

A 

tance  d  is  less  than  a  characteristic  radius  is  studied.  Specifically,  the 
intent  is  to  determine  those  time  intervals  for  which 


Rs  ±  (RrRo)2  +  4RoRi  cos2  \  sln2 


V 


2  1  2 
+  sin  —  cos 


V  + 


(IV-l) 


which  is  rewritten  as 


2  1  .2 

cos  —  sin 


V,  ,21  2  V  +  n 

— 2~  +  sin  ~2  cos  - 2 - 


(IV-2) 


where 


Rs2  -  (R1-R0>2 

4R0R1 


(IV-3) 


Unfortunately,  there  can  be  many  time  intervals  within  T^,  which  satisfy 
inequality  (IV-2),  so  that  an  accurate  analytic  solution  is  very  difficult. 
Figure  2  illustrates  this  problem  with  a  plot  of  the  right-hand  side  of  (IV-2) 
against  the  normalized  time  t/T^  for  1  =  0,  tt / 2 ,  and  Tr  radians.  The  numeri¬ 
cally  determined  percentage  of  time  that  inequality  (IV-2)  is  satisfied 
(hereafter  denoted  by  f)  is  plotted  in  Fig.  3  as  a  function  of  I  with  R^/R^  as 
a  parameter  for  that  value  of  r)  which  minimizes  f.  As  the  reader  may  expect, 
I  =  tt/2  appears  to  give  the  lowest  f  that  the  two  spacecraft  are  within  the 
specified  range.  Figure  4  shows  f  for  the  case  I  -  tt/2,  as  a  function  of 
R^/R^  for  that  value  of  n  minimizing  f.  The  irregularities  in  these  curves 
appear  to  be  caused  by  the  complicated  relationship  among  u)^,  u)^,  and  a. 

Simple  analytic  upper  bounds  are  obtained  by  considering  each  of  the 
two  products  on  the  right-hand  side  of  (IV-2)  independently.  Define  the 
times  T^  and  T^  by  the  relations: 

* 


H  was  tested  in  steps  of  tt / 18  radians. 
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Fig.  3.  Fraction  of  time  two  satellites  are  within  the  squared 
distance  a  of  one  another,  for  worse  case  phase  angle  n. 


(IV-4) 


a 


2  I  2  Vl 

cos  *2  sin  “2 — 


2  I  2  a)ZT2 

sin  2"  sin  ~2 — 


(IV-5) 


2  I  2  I 

valid  for  a  <  cos  ^  and  a  <  s*n  J*  respectively.  is  the  time  duration 

2  I  2  U)AT 

over  which  cos  sin  — rises  from  zero  to  a,  while  T2  is  the  time  duration 


2  I 


2  oj£x  +  n 


over  which  sin  cos  - 2 -  rises  from  zero  to  a.  Then  fT^,  the  time  within 


that  (IV-2)  is  satisfied,  is  upper  bounded  by 


Wr 


fTA  '  <2Tr  2  ST  V  •  * 
A 


(IV-6) 


where  the  factor  in  the  second  inequality  reflects  the  number  of  zeros 

in  an  interval  of  length  T^.  Using  the  series  expansion  of  sin  X  to  upper  bound 
and  T2  as  given  in  (IV-4)  and  (IV-5),  and  substituting  into  (IV-6), we  obtain: 


2  a 


1/2 


7 T  cos  1/2  *  7T  sin  1/2 


(IV-7) 


valid  for  all  I. 

2  I 

A  much  tighter  bound  is  easily  derived  for  the  case  a  <<  cos 

2  I  1 

a  <<  sin  -j,  by  counting  only  the  contribution  of  the  zeros  in  the 

2  I  2  U  2X2 

sin  ~2  cos  - 2 -  term  which  occur  in  the  region  where  the  cos  sin  — 

term  is  small.  Then, 


* 


The  notation  A  <  (B,C)  denotes  A  <  B  and  A  <  C. 


IV-5 


CIV-8) 


£T4  <  2 


l 


After  substitution,  we  have: 


f 


a 

sin  I 


(IV-9) 


2  i  2  I 

valid  only  for  a  <<  cos  -j*  sin  y*  This  bound  should  be  quite  tight  in  the 
region  of  small  a  and  I  close  to  ir/2.  See  Fig.  5. 


IV-6 


V.  APPLICATION  TO  SATELLITE-BORNE  SEARCH  SYSTEMS 

Consider  a  search  satellite  with  circular  orbit  radius  R^,  and  a  sensor 

with  a  detection  range  Rg  which  scans  at  random  within  a  spherical  cap  whose 

center  lies  on  the  earth-satellite  line  from  satellite  0.  The  boundary  of  the 

cap  is  determined  by  the  intersection  of  the  sphere  of  radius  R^  (R^  >  Rq)  which 

characterizes  the  satellite  being  searched  for  with  a  sphere  of  radius  R  cen- 

s 

tered  on  the  search  satellite.  By  symmetry  the  following  analysis  also  applies 
for  a  search  system  with  the  searcher  at  altitude  R^  searching  for  satellites  at 
altitude  R^,  with  a  search  system  range  of  Rg.  From  the  previous  section  it  is 
clear  that  a  satellite  being  searched  for  will  lie  within  the  cap  a  fraction  f 
of  the  total  time  T^.  In  this  section  we  will  relate  the  solid  angle  search 
rate  of  the  sensor  system  on  board  the  searching  spacecraft,  to  the  total  aggre¬ 
gate  time  required  to  achieve  the  specified  probability  of  detection  against 
the  intended  targets. 

The  sensor  system  is  characterized  by  the  angular  resolution  cell  area 

2 

0  ,  the  characteristic  integration  time  T,  and  the  single  look  probability  of 

detection  p,  which  occurs  when  the  intended  victim  satellite  is  within  one 
a 

or  more  resolution  cells  at  range  no  greater  than  Rg.  N  is  the  number  of 
parallel  detectors  operating  at  the  same  time. 

Within  the  characteristic  time  T^,  there  will  be  fT^/l  independent  oppor¬ 
tunities  for  the  sensor  to  detect  the  satellite  being  searched  for.  Within 
the  total  time  T^,  the  aggregate  cumulative  probability  of  missing  the  target 
is  given  by 

2  fTA/T 

<  <1  -  ^2  Ph>  <V-D 


Thus  for  a  total  search  time  Tg  substantially  greater  than  T^  we  obtain 


„  /.  N0  x 

Pm  1  (1  "  ~2  Pd> 

Trip 


<fTA/T)(T8/TA) 


(V-2) 


V-l 


given  by 


In  both  of  these  expressions,  the  angular  area  of  the  cap  is 

2 

steradians,  where  is  the  included  half  angle  given  by  cos  ip 


2,2  „  2 
*0  "l  -  R8 

“A 


•  It  is  easy  to  show,  using  (IV-3),  that  cos  ^  *  1  -  2a. 


Hence  for  small  a,  =  4to. 

For  the  case  where  the  detector  tests  only  a  small  part  of  the  sky  at  a 
2  2 

time,  N0  h up  is  small.  Use  of  the  exponential  approximation  yields 


-fT  2 
s  N0 


Pm  ~  e 
m 


4  to  Pd 


(V-3) 


Solving  to  obtain  the  expression  for  search  time: 


T  =  £n  (1/Pm)  (V-4) 

s  - - - 

fNft 

where  we  have  used  the  characteristic  solid  angle  search  rate: 

ft  =  e2/x  (v-5) 

In  practice,  the  probability  of  detection  for  a  sensor  varies  as  a  very 
strong  power  of  the  range  to  the  target.  In  effect,  there  will  be  a  charac¬ 
teristic  range  inside  which  the  probability  of  detection  will  be  essentially 
unity.  Outside  this  range  the  probability  of  detection  will  be  virtually 
zero.  Therefore,  for  engineering  purposes,  we  may  approximate: 


T 

s 


4to  in  (1/p  ) 
_ ro 

f  N  ft 


(V-6) 


V-2 


At 


where  f  is  computed  on  the  basis  of  the  characteristic  range  Rg  corresponding 
to  a  detection  probability  say,  .9. 

2  2 

On  the  other  hand,  if  N0  _>  thJj  ,  the  sensor  continuously  monitors  the 
cap*  Then 

fT  /T 

p  -  (1  -  p  )  8  (V-7) 

m  a 

which  gives 


£n  p 

T  -  I  _ 12 _ 

S  f  *n(l  -  PJ 


(V-8) 


with  the  restriction  Tg  _>  T^.  Generally  such  a  system  would  be  operated  to 
achieve  Tg  *  with  the  given  sensor  characteristics. 


VI .  RANGE  RATE 

In  this  section  the  range  rate  is  derived.  Equation  (III-4)  may  be 
differentiated  to  give: 


rori  r 

d(T)  |“a 


2  1  2  I 

sin  cos  *2^0)^  sin  (co^T  +  n)  sin 


(VI-1) 


Since  |d(l)  |  _>  -  R^,  we  may  upper  bound  the  maximum  range  rate  as 


•  ^0^1  r 

d<T>  -  ft  [ 


2  I  2 

0)^  cos  y  +  sin 


(VI -2) 


An  equivalent  form  is 


d(r)  < 


R0R1  (  2tt 


R0_R1  \  T0 


~  cos  I 
T1 


(VI -3) 


As  it  should,  this  expression  illustrates  the  minimum  peak  range  rate  for 
coplanar  satellites  moving  In  the  same  direction  (I  =  0°) ,  and  maximum  peak 
range  rate  for  coplanar  satellites  moving  in  the  opposite  direction  (I  =  180°) 
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